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Abstract. Stochastic mass transport models are usually described by specifying 
hopping rates of particles between sites of a given lattice, and the goal is to predict the 
existence and properties of the steady state. Here we ask the reverse question: given 
a stationary state that factorizes over links (pairs of sites) of an arbitrary connected 
graph, what are possible hopping rates that converge to this state? We define a 
class of hopping functions which lead to the same steady state and guarantee current 
conservation but may differ by the induced current strength. For the special case of 
anisotropic hopping in two dimensions we discuss some aspects of the phase structure. 
We also show how this case can be traced back to an effective zero-range process in 
one dimension which is solvable for a large class of hopping functions. 
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Stochastic transport of some conserved quantity, generically called "mass", has 
recently attracted much attention due to a large variety of applications ranging from 
microscopic (intracellular) to macroscopic (highway) traffic [1]. Important examples of 
technological interest are granular flow [2] and granular clustering [3] . From a theoretical 
point of view, these systems are challenging, since they are in general out-of-equilibrium 
and allow for phase transitions even in one dimension An example is spontaneous 
symmetry breaking and the phenomenon of condensation which happens above some 
critical mass density and corresponds to jams in traffic or aggregation in granular media. 
Nevertheless, simple models such as the zero-range process (ZRP) [5], a dynamical 
version of the balls-in-boxes model [B], or asymmetric simple exclusion processes [7], 
can capture some important aspects of these problems while remaining analytically 
solvable. 

The non-equilibrium models are defined by specifying the dynamics rather than the 
probability of a microstate. Usually one proposes the transition rates between states, 
and the goal is to predict the existence and the properties of a stationary state. This 
is the state where macroscopic observables remain constant, although some currents 
may flow in the system. In this paper we study the reverse problem. Given a steady 
state which assumes a form factorized over pairs of sites that correspond to the links of 
an arbitrary graph, we search for a class of transition probabilities which lead to this 
state. This approach is motivated by the fact that knowledge of the stationary states 
of non-equilibrium models generically facilitates the discussion of the phase diagram, 
because a number of observables can be calculated analytically. For instance, the ZRP, 
deflned in terms of particles hopping between sites of a lattice and interacting only if 
they are at the same node, has a steady state that factorizes over the sites of a lattice, 
or more generally, over nodes of an arbitrary graph. The factorization allows for a 
convenient mathematical treatment. A generalization of the ZRP that leads to pair- 
factorized steady states (PFSS) was proposed in [8] for a one-dimensional ring topology. 
It was shown that nearest-neighbor exponentially suppressed interactions plus some 
extra "pinning" (ZRP-like) potential result in a condensate that is spatially extended. 
In [9] the shape of the condensate was derived. It was also shown that the scaling of 
the extension of the condensate with the system size can be tuned via an appropriate 
competition between local and ultralocal hopping interactions. 

We shall show in this paper that the ZRP and the PFSS on a ring (where PFSS 
here should be understood as the corresponding processes leading to PFSS) are special 
cases of a more general setting. Beyond that, we shall consider non-local processes on an 
arbitrary graph, or processes with anisotropic hopping in two dimensions. In the latter 
case we shall show that it can be dimensionally reduced to a ZRP in one dimension with 
weights that contain the information on the pair-factorized stationary behavior in the 
second dimension. Therefore former results on PFSS on a one-dimensional ring topology 
^ can be used to derive features of the condensation transition in the anisotropic case. 

The model. We consider a connected, undirected but otherwise arbitrary graph 
with N nodes (sites), and node degrees ki, ...,kN- We place M particles of unit mass 
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on the nodes of this graph. "Particles" here stands for a generic mass that is involved 
in the transport process and has "bosonic" properties in the sense that mj > particles 
may be assigned to the same site i. The distribution of occupation numbers of nodes 
is denoted as m = {mi, mTv}- The dynamics is defined as follows. We pick up 
a randomly chosen node i, and if it is not empty, a single particle departures with 
probability Ui{mi \ . . .), where the dots stand for occupation numbers of other sites (in 
general not necessarily nearest neighbors of node i). Next, the particle chooses a target 
site j with probability Wij = W{i — > j). The transition matrix W may be arbitrary, 
with the only assumption that all Wij > and J2j ^ij = 1 ^- The hopping 

event is thus spht into two steps: the departure from a site, determined by the function 
u, and the choice of destination site, determined by the rates Wij. 

We will next derive under which conditions on the hopping rate and the transition 
matrix the system reaches a steady state that assumes a pair-factorized form 

N / N \ 

w{m) = Yl gijimi, m,) JJ k^^6 m, - M . (1) 

i=l \i=l J 

Here denotes pairs of nearest neighbors, gij{mi,mj) is a symmetric but otherwise 
arbitrary weight function that may depend on the link and the (5-function that 

ensures conservation of the overall mass M and will be dropped below. The factors 
k^* are slit off for convenience. As we shall prove next, a possible class of hopping 
rates (although not the only possible one, since, for example, the symmetry condition 
on g{m, n) was released in jSj) is given by 

n "-^'T'-?' . (2) 

where Af{i) denotes the set of neighbors of node i provided the transition matrix Wij 
satisfies the condition 

k = J2 W,^k, . (3) 
j 

Proof. Let us start with the balance equation for a conserved probability current at 
each node i, this is a condition that is necessary for having a steady state: 

Ui{mi\ . . .)w{m) = '^W{j i)uj{mj + 1| . . . , - 1, . . .) 

i 

X w{mi, . . . , m-i - 1, . . . , mj + 1, . . . , mA?) . (4) 

Equation (|4]) says that the probability with which a particle leaves a site i with m^ 
particles in a configuration given by occupation numbers m should be the same as the 
total probability that node i receives one particle from any site j that was formerly 
in a configuration with nij + 1 particles at site j and rrii — 1 particles at site i. The 
individual probabilities for a hopping event from j to i are given by W{j — > i). As a 
sufficient condition, the probability current is conserved if ^ is satisfied individually 
for any i E {1, .., A^} and any set of occupation numbers m. Now, dividing ^ by w{m), 
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using the definition ([T]) of the steady state, and inserting ^ as an assumed form of 
Ui{mi \ . . .), we cancel common factors in the ratio w{...)/w{...) on the right-hand side 
and expand it to obtain 

gij{mj,mi-l) yj gaj{mj,ma) ^ 



Ui[mi\ 



gijiruj + 1, mi - 1) ^AJ.^,^ g^.^j^. + 



X 



-i-r gbijrni - I, nib) gijim - l,mj + 1) -|-|- gcjjmj + l,mc) 



Note that the first bracket is just Uj, rewritten as a product of ratios of g, with the part 
for a = i written out exphcitly. Similarly we express the ratio of w{...) /w{...) in terms 
of the weight functions g and explicitly split off the factor for b = j. The first product 
(over a) and the last product (over c) as well as the terms gij cancel out. We end up 
with 

where the product over b now includes all nearest neighbors of i. The expression in the 
bracket is just Ui{mi \ . . .). Therefore, equation ^ is satisfied if 

k = Y, W,,k, . (7) 



This is precisely (|3j) which completes the proof. □ 

Before we come to examples, we want to add a remark on a useful possible 
reinterpretation of the involved links included in Wij. Since the stationary state ([T]) 
only refers to the occupation number distribution and depends on the hopping function 
■u, but not on the transition matrix Wij, by making different choices on Wij we can tune 
the current of particles in the system. Moreover, the set of non-zero Wij does not need 
to coincide with the set of existing links in the graph. This may be seen as if we had two 
graphs: an undirected one that specifies which pairs of nodes enter the steady state ([T]), 
and a directed and weighted graph of allowed transitions, defined by Wij = W{i — > j). 
This fact can be used in numerical simulations of the steady state. Suppose that we have 
a process which, although being out of equilibrium, leads to ([T]). If we are interested only 
in static properties of the steady state and have to turn to Monte Carlo simulations as 
the model is not fully analytically accessible, we can change the dynamics (by changing 
Wij) in order to speed up the convergence towards the stationary state, or to reduce 
the autocorrelation time. An example is to turn the local hopping of particles with 
given hopping rates into a non-local update scheme by choosing two random nodes and 
moving the particle between them, accepting or rejecting the move according to the 
rules of the Metropolis algorithm. This can be viewed as specifying Wij = 1/{N — 1) for 
all i,j, i.e, as making the "transition" graph a complete graph. One can check that such 
a choice obeys ([s]), and since the probability of moving the particle along each link is 
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the same in both directions, there is no macroscopic particle current. More generally, if 
Wij = 1/ki for every link equation ^ is fulfilled and the current vanishes, so the 

system is actually in equilibrium. This method can in particular be used to simulate the 
PFSS on a ring mentioned above. In the condensed phase, effectively non-local updates 
via the Metropolis algorithm accelerate the convergence to the stationary state, since 
it facilitates the melting of local, separated, metastable aggregates to finally merge to a 
single condensate. 

Let us now discuss some examples. Choosing gab{m,n) = p{m), i.e., making the 
two-point weight independent of the second argument and being the same for all links, 
one obtains Ui{mi\ ■ ■ ■) = u{mi) = p{mi — l)/p(mj). Here u{m) depends only on the 
occupation number at the departure node. If one now chooses Wij in any way that 
satisfies ([3]), the stationary state w{m) = p{mi) ■ ■ -pimi^) is precisely the steady state 
of a ZRP on an arbitrary network [5]. 

A more complicated example is a one-dimensional closed chain so that ki = 2 for 
all sites i, and gab{^, n) = g{m, n) being independent of the link, but depending on 
occupation numbers at its both ends. Equation ([2]) then leads to 

Ui(mi|mi_i,mi+i) = p{mi,mi+i)P{mi,mi_i), (8) 

where P{m,n) = g{m — l,n)/g{m,n). This is essentially the PFSS considered in jS]. 
The only difference is that no assumption on the symmetry of g{m, n) was made there. 
It has been shown in [S] that when 

g{m,n) = exp[-J|m - n\ + U{6ra,o + 5„,o)/2], (9) 

in the thermodynamic limit the system exhibits a condensation transition above some 
critical density of particles p = M/N which is a function of the parameters J, U. The 
nature of this condensate is different from an analogous phenomenon in the ZRP, because 
it emerges due to nearest-neighbor interactions and not due to the on-site potential p(m) 
as in the former case. As a result, the condensate is extended over ~ \/N sites, see [8] 
and ^ for details. 

Another simple special case is a square lattice with periodic boundary conditions 
on which we can specify hopping rates that factorize over adjacent sites. Let us denote 
a site on this two-dimensional torus by its coordinates From ([s]) we obtain that 

W{i - 1, J ^ z, j) + W{i + 1, J ^ z, j) 

+W{i,j + 1 z, j) + W{i,j - 1 z, j) = 1, (10) 

so that the probabilities for ending up at target site (i, j) should add up to 1. Similarly, 
the probabilities for leaving a site (z,j) as departure site should add up to 1, that is 

W{z,j ^ I + 1, j) + W{i,j -> I - 1, j) 

+W{i,j ^ 2, J + 1) + Wit, J ^ J - 1) = 1. (11) 



Equations (10) and (11) are the only conditions on transition probabilities. They have 
many solutions, which differ by the degree of current anisotropy, e.g., the current may 
flow only in horizontal, vertical, or in both directions. 
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Anisotropic interactions in two dimensions. Now we shall describe a case which 
(from a theoretical point of view) is of particular interest due to its analytic solvability 
and the existence of two phase transitions. Let us assume that we have a two-dimensional 
lattice with N = L x L sites, periodic boundary conditions and the following steady 
state: 

w{m) = Y\Y\g{mi,j,rni,j+i) S I ^rriij -Ml. (12) 

i=l j=l \ i,j / 

This implies that g{m, n) lives on links only in one, say vertical, direction. On horizontal 
links, the weight is assumed for simplicity to be constant and equal to 1, but in general 
one may choose any weight /(mj j) depending only on a single site. The hopping rate 
assumes now the form 

uinii A . . .) = ^^^^ ' M i^ jAM ' 13 

g{mij,mij_i) g{mij,mij+i) 

which is very similar to (p|, but the particles can now jump either in horizontal or vertical 



direction. The current is defined by any W obeying equations (10) and (11). Although 
we will not address dynamical issues here, for definiteness let us assume that the particles 
can jump only to the right with probability W{i,j i,j + 1) = p or to the top with 
probability W{i,j — > z + 1, j) = I —p, while W{i,j — > i, j — 1) = W{i,j — > i — 1, j) = 0. 
Moreover, let the weight g{m,n) be chosen according to (|9]). Based on our knowledge 
on the separate ZRP and PFSS processes in one dimension, we would naively expect the 
following scenarios to be true: a) below a certain mass density, particles distribute over 
the two-dimensional grid homogeneously like a liquid; b) above some critical density, 
as a remnant of the interaction-dependent hopping in vertical direction, an extended 
condensate forms along each vertical line; c) above some critical density, as a remnant of 
the ZRP in one dimension, the overall condensate gets localized along a single vertical 
line, along which the form of the condensate is determined by the one-dimensional 
PFSS; d) the condensate stays localized along the horizontal axis, but looks like a liquid 
along the condensate-carrying vertical line. Figure [T] shows snapshots of Monte Carlo 
simulations for J = U = 1 and increasing density p. For our choice of weight functions we 
see that scenarios (a), (c) are realized below and above some pi, respectively. Scenario 
(d) happens for some p2 > pi, while scenario (b) is excluded. We shall derive next 
why only these scenarios can happen. We will argue that pi is the critical density of 
the corresponding PFSS in one dimension and also calculate p2- Note that the whole 
system does not simply split into a direct product of a ZRP in one dimension and a 
PFSS in the other dimension, since both subsystems are coupled via the overall mass 
conservation. 

First we observe the following. In the steady state we can treat the probability 
w(m) as the weight of a microstate of a system being in equilibrium, therefore we can 
formally write the partition function of the system in the canonical formulation as 

L L / 

z,2d{N =L\M)= XI n n ^i^^'i^ "^^i - ^ 

{TTiij} i=l j=l \ i,j 
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Figure 1. Snapshots of Monte Carlo simulations of the anisotropic system with 
L = 200, J — U — 1, for various densities p = 0.2 (upper-left), 0.3 (upper-right), 
0.35 (bottom- left) and 0.5 (bottom-right). For p = 0.2 the system is in the liquid 
phase. For p — 0.3 it is slightly above pi « 0.24 — a condensate emerges. For 
p = 0.35 the condensate is fully developed. For p — 0.5 > p2 ^ 0.45, the borders of 
the condensate merge and the whole vertical line is uniformly covered. 



M ML 
Mi=0 Ml=0 i=l 

where 

M ML 

Z,UL, M)=J2 Y.]1 ( M - J] ) (15) 

mi=0 mL=0 i=l 



is the partition function of a Id PFSS. Since (14) has the same functional form as the 
partition function for the zero-range process |TD] , the partition function Z^id may be seen 
as the weight p{m) = Zcid{L,m) that is now associated with the total mass m = Mi 
along the ith vertical line. 

In [9] we have shown that for m > pcL, where pc is the critical density for 
condensation in the one-dimensional system, Zcid{L,m) behaves as ~ exp{—c^/rn) with 
some c > 0. This means that for large m, the hopping rate of the corresponding ZRP, 
u{m) = p{m — l)/p{m) = ZddiL, m — l)/Zcid{L, m), behaves as u{m) = 1 + c/ (2y/rn). 
For such a hopping rate it is known [TU] that the ZRP exhibits a condensation 
transition. The condensate occupies a single site and the fluid-phase distribution is 
a stretched exponential distribution. Since m denotes now the mass along a vertical 
line, the mapping from the ZRP back to the anisotropic model allows us to predict the 
spontaneous symmetry breaking into the state with a condensate so that one of the 
masses Mj will grow to pick up all the difference AM = M — piL"^ = L^(p — pi) between 
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the total mass M and the mass piL^ in the critical background. So, differently from 
scenario (b), not each line carries its own condensate, but the condensate is localized 
onto a single line if there is a condensate at all, while all other lines have a mass density 



close to the critical value. To predict the critical density pi we observe that ( 14 ) in the 
grand-canonical formulation, 

L 

Zg,,{L,z) = J2Zc2d{L,M)z'" = J2llZ,UL,M,)z^'' 

M {Mi} i=l 

^Z,irf(L,M)z^^j , (16) 

becomes just a power of a grand-canonical partition function in one dimension. Since the 
condensation transition is determined by the non-analytic behavior of Zg2di the critical 
density pi must be the same as Pc of the one-dimensional PFSS. For J = U = 1 it reads 
pi ^ 0.24, as follows from [9]. 

In order to check how the particles are distributed along the line which carries the 
condensate, in particular to determine the shape and extension of the condensate, we will 
again make use of results of ^ for the one-dimensional system. It is proved there that 
for sufficiently large systems, the condensate has a quasi-parabolic shape with very sharp 
borders. Its width W grows proportionally to VAM, so that here W = woLy/p — pi 
with some constant Wq that depends only on J, U from ^ and reads Wq ~ 2.2 for 
J = U = 1. This is precisely the shape seen in figure [T} for pi < p < p2- When 
the density exceeds p2 = I/wq + Pi ~ 0.45, the width W becomes equal to the linear 
size L. Since fiuctuations of occupation numbers in the condensate can be neglected 
in the thermodynamic limit for large systems and p > p2 there are no empty sites. 
The ultralocal weight e^^^"^'°~^^"'°^^'^ becomes equal to one and hence g{m,n) effectively 
behaves as exp{—J\m — n\) as if U were zero. We may therefore see the line carrying the 
condensate as a new PFSS with U = 0. From [8] we know that such a system is always in 
the liquid state. This means that at p = p2 the system undergoes a second (geometric) 
phase transition to a state, in which both borders of the condensate merge and the 
particle distribution along the condensate line looks uniform apart from fluctuations, 
see flgure[T| which corresponds to scenario (d). We call it a geometric phase transition 
as the condensate percolates all over the line. It is not just a flnite-size effect of merging 
borders, because p2 remains flnite also in the thermodynamic limit L ^ oo. Thus the 
condensate changes its shape from quasi-parabolic in vertical and localized in horizontal 
direction, to a homogeneous distribution vertically, but still localized horizontally. Any 
transient remnant of the former condensate peak gets rapidly washed out towards a 
homogeneous distribution. 

In summary, we proposed a class of hopping rates that leads to pair-factorized 
steady states on arbitrary graphs. The proof holds in particular for non-local and 
inhomogeneous hopping rates. As an example we studied an anisotropic two-dimensional 
system with PFSS, for which we predicted the onset of condensation including the shape 
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and the scaling of the condensate by dimensional reduction of the two-dimensional 
system to an effective zero-range process in one dimension. 
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